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We describe a system for the compensation of time-dependent stray magnetic fields using a 
dual channel scalar magnetometer based on non-linear Faraday rotation in synchronously opti- 
cally pumped Cs vapour. We detail the active control strategy, with an emphasis on the electronic 
circuitry, based on a simple phase-locked-loop integrated circuit. The performance and limits of the 
system developed are tested and discussed. The system was applied to significantly improve the 
detection of free induction decay signals from protons of remotely magnetized water precessing in 
an ultra-low magnetic field. 
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I. INTRODUCTION 

High sensitivity magnetometry has a wide range of 
applications in both fundamental and applied research. 
Nowadays, the best performance in terms of sensitivity is 
provided by Superconducting Quantum Interference De- 
vices (SQUIDs) and by Atomic Optical Magnetometers 
(AOMs) operating in the Spin-Exchange Relaxation-Free 
(SERF) regime. In both cases, the impressive limit of 
1 fT/VHz has been surpassed experimentally, when op- 
erating in volumes carefully shielded against the environ- 
mental magnetic noise. [l| 

Scalar magnetometers operating in non- vanishing mag- 
netic helds have also been reported to have a sensitiv- 
ity as good as a few tens of fT/VHz when operating 
in shielded environments. This value degrades to about 
1 pT/yTIz in unshielded volumes. 

The precise determination and control of magnetic 
fields is a key factor in atomic physics experiments [2|-|j] , 
fundamental research (such as determination of the neu- 
tron electric dipole moment [El, HI), ultra-low- field (ULF) 
NMR and imaging Ml , and the development of biomag- 
netic diagnostics |1,E|. 

Field control is commonly attempted through both 
passive shielding and active compensation approaches, 
and many solutions are reported in the literature. Differ- 
ent applications require that the shielding/compensating 
systems have different levels of noise rejection and re- 
spond in different frequency ranges. 

Passive shielding can be achieved by surrounding the 
measurement volume with soft-iron/large permeability 
material, or with thick conducting material (or super- 
conducting layers) in which eddy currents are induced. 
The superconducting approach has been studied and de- 
veloped since the 1970s |9h11| and yields the highest field 
stability. 
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The efficiency of passive shielding decreases at low fre- 
quency, thus shields can be profitably coupled with ac- 
tive compensation systems to improve the performance 
in the Hz and sub-Hz range. Active compensation has 
also been proposed and studied for use in compensat- 
ing peaked spectral components of the stray field, such 
as those at 50 Hz or 60 Hz due to power lines. Active 
compensation of deterministic disturbances has also been 
attempted with a feed- forward approach [12]. Neverthe- 
less the most recent literature refers to feedback systems, 
which arc suitable for attenuating both deterministic and 
random components. 

Active compensation with a suitably extended band- 
width can also be applied in stand-alone configurations 
(i.e. not combined with passive shielding). This is of 
interest, for instance, where easy access (no geometrical 
constraints) is required, and/or when large magnetically 
clean volumes are needed, which would make the passive 
shield bulky and very expensive. 

Naturally, the specifications (bandwidth and attenua- 
tion factor) and the final performance of such active set- 
ups depends on the feedback design and on the sensors 
used to generate the error signals. In the latter respect, 
a large variety of sensors were reported in the literature, 
ranging from the mechanical ones (a sort of compass!) of 
decades ago [l3[ to Hall sensors fl3 . magnetoresistances 
IJH, fluxgates [H3, SQUIDs~M[II, SERF- AOMs 
19| , NMR based magnetometers [2fJ . 

A further approach frequently used to counteract mag- 
netic noise relies on the fact that noise sources, being at a 
large distance from the measurement volume, produce a 
rather homogeneous disturbance. Thus the use of differ- 
ential (gradiometric) sensors make it possible to reduce 
noise when measuring signals produced by nearby sources 
[2lj . so that the common- mode signal can be attributed 
to magnetic noise and rejected. 

In this paper we present the results obtained with a 
dual channel scalar Cs AOM operating in /iT field, which 
was used to stabilize the modulus of the bias field. This 
set-up makes it possible to detect the low-level time- 
dependent anomalies generated by small sources located 
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close to one of the two sensors composing the differential 
magnetometer. Both sensors work in the same bias mag- 
netic field, which is also the field in which the anomaly 
source is immersed. The stability of the bias field plays a 
key role when the magnetometer is used to detect NMR 
free induction decay: stabilizing the field makes the nu- 
clear precession occur at a more stable frequency, thus 
facilitating its analysis and long lasting averaging. 

While the characteristics of the dual channel magne- 
tometer and of its applications to ultra-low-field NMR 
detection [22|, [23| and magneto-cardio signal detection 
pij are reported elsewhere, here we discuss the perfor- 
mance of a bias field stabilization system, focusing in 
particular on the effects of residual magnetic noise in our 
simplified configuration, in which a pair of scalar sensors 
are used to counteract the noise in the field modulus by 
acting on a single component. 

In the case of a self-oscillating magnetometer, con- 
sidered as a field-to-frequency converter, when the field 
compensation is controlled by a voltage, the entire setup 
acts as a voltage-controlled oscillator: the whole com- 
pensation system can be regarded as belonging to a well 
known and extensively studied class of electronic circuits, 
phase- locked loops (PLLs). From a practical and tech- 
nical point of view, this analogy has two advantages: a 
complete feedback theory related to the PLL has already 
been developed and assessed, and the fact that a variety 
of integrated circuits are commercially available simpli- 
fies the task of building high quality phase comparator 
circuits of several kinds, and designing loop filters. The 
system described in this paper uses a commercial PLL 
integrated circuit and a micro-controller chip (PIC) used 
as a programmable divider to generate an adjustable ref- 
erence frequency. We foresee that further improvements 
could be achieved by using digitally controlled PLL de- 
vices, which contain different kinds of phase comparators 
and programmable loop filters [251 ] . 



II. EXPERIMENTAL APPARATUS 

The experimental apparatus is sketched in Figs. [1] and 
[2j Two sealed cells containing Cs and 90 Torr of Ne as a 
buffer gas are heated to about 45° C. The two cells are 
displaced by 7 cm along the x direction. Each cell is illu- 
minated by two overlapping laser beams, both attenuated 
down to a few //W/cm 2 . One of them (the pump beam) 
is circularly polarized and broadly frequency modulated 
so to make it pass periodically in and out of the Doppler- 
broadened atomic resonance. The other beam (the probe 
beam) is linearly polarized and unmodulated. The bias 
magnetic field (Bqz) is perpendicular to the beams' prop- 
agation axis (y). Provided that the pump beam excites 
atoms synchronously with their Larmor precession, a pro- 
cessing macroscopic polarization of atoms is induced in 
both the cells. The probe beam polarization undergoes 
a time-dependent Faraday rotation of the polarization, 
which is detected by means of two balanced polarime- 



ters. The modulation frequency of the pump laser can 
be scanned to detect and characterize the resonance, or 
be produced by a pulse generator triggered by the po- 
larimetric signal of the main channel, as represented in 
Fig. [1] The latter configuration makes the system a self- 
oscillating magnetometer. The plots show the resonance 
shape detected in scanned mode and the dependence of 
the oscillating frequency as a function of the loop de- 
phasing introduced as a feedback delay. The merit factor 
of the resonance (Q = vq/ Avfwhm) is about 130 and 
is consistent with the maximum slope of the lower plot 
(7rrad/110Hz). This value was optimized by adjusting 
several experimental parameters such as cell tempera- 
ture, laser intensities, and probe detuning. The value 
achieved is consistent with other optimized results re- 
ported in the literature [26]. Additional details of the 
self-oscillating setup are available in [22|, [23| . 

A. Coil dc supplies 

Fig. [5] represents (not to scale) the relative positions of 
the atomic sensors and compensation coils. A set of three 
mutually orthogonal, large (180 cm sides) Hclmholtz 
pairs are used to control the magnetic field. Two pairs 
are used to cancel out the local magnetic field in the hor- 
izontal plane (a; and y pairs) , while the third pair is used 
to partially compensate the z component, thus setting 
the value of the bias field to about 4 /iT. Three separate 
current generators are used to supply the compensation 
coils. The current supply for the z pair is analogically 
controlled by high quality potentiometers, while the cur- 
rents for the x and y pairs are digitally (and automati- 
cally) controlled. The procedure used to zero the trans- 
verse components of the field is based on the minimiza- 
tion of the self-oscillating frequency: the PLL circuit is 
deactivated and a simplex procedure adjusts the currents 
in the x and y coil pairs while evaluating and minimizing 
the atomic Larmor frequency. 

B. Active phase-locked loop compensation 

The analogue signal generated by one balanced po- 
larimeter is preamplified, filtered and compared to a 
threshold value, in order to produce a synchronous square 
wave which triggers (possibly after frequency scaling) the 
self-oscillating magnetometer. The same signal is phase- 
compared with a stable reference signal in order to pro- 
duce the error signal for the PLL. The phase compari- 
son is performed by phase-comparator- 1 (PhC-1, exclu- 
sive or) or by phase-comparator-2 (PhC-2, edge-triggered 
three-state circuit) of a CD4046 integrated circuit. 

A chain made of a one-pole/one-zero active filter fol- 
lowed by an adjustable-gain linear amplifier and by a 
voltage-to-current converter, closes the loop, sinking a 
small part of the current supplying the z Helmholtz pair, 
see Fig. [3] 
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FIG. 2: Schematic of the apparatus (not to scale). The trans- 
verse field components, B x<y , are compensated by Helmholtz 
coils (HC) supplied by numerically controlled current gener- 
ators. The B z component (bias field) is analogically adjusted 
and actively stabilized. The symbols • and o represent the 
cells of the main arm and secondary arm, respectively, while 
represents the sample. For the sake of clarity, one of the 
transverse HCs, and the gradient compensation quadrupoles 
are not shown. 



activated the compensation system: the 50 Hz sidebands 
are reduced to the broadband noise level (i.e. by about 
60 dB) in the case of the main arm and by 42 dB for 
the secondary arm. In addition, for both the arms the 
carrier peaks emerge from much weaker pedestals, indi- 
cating a significant rejection of low frequency noise. A 
quantitative analysis of such noise reduction is made in 
the following. 



FIG. 1: (a) Schematic of the self-oscillating optical atomic 
magnetometer (main arm). A balanced polarimeter detects 
the atomic precession, generating a signal that triggers the 
synchronous optical pumping. The same signal is phase- 
compared with a stable reference signal to generate the er- 
ror signal, which is fed back to stabilize the field (see Fig.[3|. 
The (b) plots show the in-phase and quadrature signals ob- 
tained in scanned mode. The (c) plot shows the relation- 
ship between loop dephasing and oscillation frequency in self- 
oscillating mode. 



III. PERFORMANCE ANALYSIS 

The self-oscillating magnetometer converts the modu- 
lus of the magnetic field to frequency, so that any modula- 
tion in the modulus appears as a pair of sidebands in the 
power spectral density of the magnetometric signal. The 
left plots in Fig.@] show typical spectra recorded without 
any shielding or active compensation. Strong sidebands 
appear at the power line frequency and its multiples, as 
a result of significant magnetic noise produced by trans- 
formers and other electrical devices. The amplitude of 
such sidebands varies in time, and in the case of the spec- 
tra shown amounts to -10 and -13dB c for the main and 
secondary arm, respectively. The right plots in Fig.0] 
were recorded in the same conditions, apart from having 



A. Short-term 

The analysis of the short term performance consists 
in the characterization of the residual phase noise in the 
polarimetric signals. The relatively low-frequency oper- 
ation of the magnetometer makes the direct digitization 
of these signals possible, as well as the implementation 
of a full numerical procedure for phase noise evaluation. 

The data analysis reproduces, using numerical tools, 
the typical setup of double balanced mixers (DBM), 
which are commonly used in the phase noise characteri- 
zation of RF and microwave devices (see [27l . [28j | for an 
overview and more in depth discussion). Pair of sam- 
ples are digitized at 50kS/s (or less) with an amplitude 
resolution of 16 bit. Each sample is a trace containing 
2 13 = 8192 data points. In these traces two out of the 
three signals available are recorded: the main arm polari- 
metric signal, the secondary arm polarimetric signal, and 
the PLL reference. Consequently, noise characterization 
can be achieved for the relative phase of the two arms or 
for the absolute phase of each individual arm. 

Once the digitized traces are transferred to the com- 
puter, the main single tones of the two traces are evalu- 
ated by means of a Lanczos-Grandke procedure [1^, [3(| , 
in order to extract their frequency and relative phase and 
to normalize their amplitude. Both traces are then in- 
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FIG. 3: Schematic of the loop filter and gain adjustment cir- 
cuit. The phase comparator (PhC) signal is filtered by OA1, 
amplified (and optionally inverted) by OA2, and converted 
to current by OA3 and Ql. OA4 supplies a local reference 
biased 3V above the ground. Two green LEDs are used as a 
monitor for the Ql work-point. The Ql collector sinks part 
of the Helmholtz coil current, when the closed-loop operation 
is selected. 



terpolated by means of FFT/zero-padding/FFT^ 1 tech- 
nique, up to an apparent sampling rate of about lMS/s 
or more. The two digitized signals are then relative phase 
shifted up to a dephasing of 7r/4. The precision of the 
phase shifter is limited by the sampling rate 1 / dt' of the 
interpolated trace (<5</> = l-Kv^dt'}. This amount sets an 
upper limit to the accuracy of the phase measurement de- 
scribed below. The resulting quadrature traces are then 
multiplied element-by-element. The procedure described 
so far implements the phase-shifter and the mixing cir- 
cuit of the DBM. The low pass filter which follows the 
mixer in DBM setups is replaced here by a single tone 
identification and subtraction routine, tuned to the sec- 
ond harmonic of the previously detected frequency The 
residual signal can thus be modeled as 



g(t) = sin[27rz/ £ + 0i (t) + 5(f)] cos[27w i + 4>2(t)] 
=-(sin[2 x 2nv t + 0i(t) + <f> 2 (t) + 5<t>] 
+ sin[<ty + 1 (t)-0 2 (t)]) 



(1) 



which, after having filtered out (actually subtracted) the 
component at 2i/q, can be approximated as 
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FIG. 4: Power spectrum of the polarimetric signal. In (a) 
and (b) no compensation is used and 50 Hz noise sidebands 
are clearly visible, with an amplitude of about —10-: — 13dB c . 
The plots (c) and (d) are obtained with active compensation. 
Plots (a) and (c) refer to the main arm, while (b) and (d) refer 
to the secondary one. Plots (a) and (c) show that the 50 Hz 
sidebands disappear, meaning that an extinction of 60 dB is 
attained in the main arm. In the secondary arm (plots (b) 
and (d)) the 50 Hz sidebands are attenuated by 42 dB. Thanks 
to the rejection of low frequency noise, in traces (c) and (d) 
the carrier peak emerges from a much weaker pedestal with 
respect to (a) and (b). 



9 ^ ' ~ 2 ~ 2 2 

.(2) 

where the two approximations consist in the lineariza- 
tion of sin(A0) in the hypothesis A0 = <f>\ — 02 "C 1, 
and in the assumption that 5<f> — 0, respectively. The 
phase difference A(j>(t) is finally analyzed by repeated 
FFT evaluations of its power spectral density, which are 
then averaged to highlight the main features of the spec- 
trum. The above mentioned quadrature error S<p implies 
a systematic error (1 — cosScf)) 2 ~ <50 4 /4 on the spec- 
trum scale, which is negligible provided that the inter- 
polation produces a good level of oversampling, e.g. we 
used dt' 1/(10(H) 60mrad. 

The phase noise spectra recorded by comparing differ- 
ent couples of traces are plotted in Fig. [5] Absolute phase 
noise in the main arm (o) appears together with the rela- 
tive phase noise of the two arms, recorded in locked-ficld 
(A) and unlocked-field (•) conditions (right axis units). 
The relative phase plots are used to infer the differential 
magnetometric sensitivity, which is evaluated by convert- 
ing the power spectral density of the phase noise to mag- 
netic field units (left axis units): 
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dip v iQ 
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(3) 



where 7 ~ 21.97 x 10 9 (rad/s)/T is the gyromagnetic 
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FIG. 5: Residual absolute phase noise in the main arm (o), 
and relative phase noise between the two arms in locked (A) 
and unlocked (•) field conditions (see text for details). The 
right and left vertical axes refer to phase noise and magnetic 
sensitivity, respectively. The conversion between them follows 
eq.rj 



by the substantial agreement between the plots (•) and 
(A) in Fig. [5] In short, the common mode noise repre- 
sented by magnetic fluctuations is counteracted by the 
system, while phase fluctuations originating from other 
noise sources are compensated by magnetic field adjust- 
ments in the main arm. In this way, part of the noise 
detected in the main arm is transferred to the secondary 
arm through an additional magnetic noise. 

The data shown in Fig. [5] and in Fig. [4] were obtained 
using as a phase-comparator the PhC-1 (exclusive or). 
We observed that, in spite of a larger residual noise at 
2vq in the error signal (partially transferred to the coil 
current), PhC-1 gave a lower background in the phase 
noise than PhC-2. This feature is consistent with the well 
known greater sensitivity to noise of PhC-2. A cleaner 
polarimetric signal would probably result in favourable 
conditions for the use of a PhC-2 instead. 



B. Long term 



factor of Cs. 

The results shown in Fig. [5] demonstrate an ultimate 
limit of the residual noise of less than 1 pT/VSz in the 
main arm, and of 2.5pT/\/Hz in the secondary one. 
Concerning the noise at 50 Hz, residuals of 4 pT/VHz 
and lOpT/vTIz appear in the two arms, respectively. 
This evaluation is consistent with the rougher estimate 
derived from Fig. [4] It is important to point out that 
the 50 Hz noise varies in time, thus preventing accurate 
comparisons. The peak at 150 Hz is not a third-order 
sideband, but rather a consequence of the 50 Hz noise 
anharmonicity, and is only partially attenuated in both 
arms as it is located close to the cut-off frequency of the 
loop filter. A spurious peak at 100 Hz, is also occasionally 
observable in the phase noise of both arms. This cannot 
be explained in terms of a second order sideband (as it is 
definitely larger than estimate Jf (M)/Jf(M) rj M 2 /16 
at the modulation index M inferred from the first side- 
band amplitude). This is, for example, the case of the 

Fig.ru (d). 

The phase noise spectra can be compared to the noise 
level expected from the photodetectors and photocur- 
rent amplifiers. The photocurrent shot noise, tran- 
simpedance Johnson noise, and preamplifier input cur- 
rent noise all contribute with similar amounts, resulting 
in about 135^V/-\/Hz at the output of each arm. This es- 
timate, with an in-phase signal V x ~1 V, leads to a white 
phase noise as large as Sip 2 « SV 2 /V 2 « 20nrad 2 /Hz, 
consistently with the values that appear in Fig. [5] 

The stabilization system reduces such phase fluctua- 
tions by a factor of 4 in the main arm. This is achieved 
by the desired and appropriate reaction to the magnetic 
noise and to an unwanted reaction to other noise sources. 
The latter is at the expense of an increase in the mag- 
netic noise in the secondary arm. It is worth noting 
that no additional noise is introduced, as demonstrated 



Even in ideal conditions, a test of the long term sta- 
bility of the compensation system described, cannot be 
performed by analyzing the drift or the Allan variance 
of the self-oscillating frequency: such a test would in- 
stead characterize the relative stability of the reference 
oscillator and the data acquisition clock. 

Evaluation of the long term stability, i.e. the amount 
of the noise contribution in the sub- Hz sub-mHz range, 
requires knowledge of the extent to which a certain value 
of the (locked) oscillation frequency corresponds to a 
given value of the magnetic field. In fact, while the 
self-oscillation frequency is ideally set by the Larmor fre- 
quency of atomic spins and thus by the modulus of the 
magnetic field, several subtle effects may induce slight 
deviations. These effects have been studied intensively 
as they are commonly considered as sources of system- 
atic errors in precise magnetometry. Should they vary 
in time due to drifts in the atom-light interaction con- 
ditions, an apparent field drift would appear and would 
be inopportunely compensated by the active stabilization 
system. In conclusion, the long term stability would be 
set by the drifts of the systematic errors. 

We point out several potential causes of drifts of sys- 
tematic errors, among which the most relevant are the 
cell temperature and the laser junction currents and tem- 
peratures, with their well known effects on laser tun- 
ing. In particular, drifts in pump laser detuning result 
in slight phase deviation of the synchronous pumping, 
thus producing drifts in the self-oscillating frequency, in 
accordance with the behaviour shown in Fig.[T] plot (c). 

The relationship between loop dephasing and oscilla- 
tion frequency is a delicate feature of self-oscillating mag- 
netometers [3l[ . The loop dephasing is set by a time de- 
lay in the loop, whose conversion to dephasing depends 
on the operating frequency. For this reason, stabilizing 
the frequency makes the magnetometer more reliable and 
robust. 
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It is worth noting that the long term drifts limit the 
performance of the magnetometer to a greater or lesser 
extent depending on the specific application. We con- 
sidered applications in magnetocardiography [23j and 
in ULF-NMR remote detection |22j . In the first case 
only drifts occurring during a cardiac epoch are relevant, 
which are definitely negligible compared to the signal dy- 
namics. In the case of NMR detection, the low signal- 
to-noise ratio makes long lasting averages necessary, so 
appropriate measurement requires that the field drifts oc- 
curring during the whole average process causes the nu- 
clear precession frequency to drift negligibly compared 
to the NMR resonance width. We considered averaging 
intervals ranging from several minutes to several hours, 
and the nuclear resonance detected (at v n s» 180 Hz) had 
a width of Av n = 1/(2ttT 2 ) wlOOmHz. The stability 
required can thus be written as AB/B — Av n jv n <C 
5 x 10 , to prevent an apparent reduction of T2 due to 
inappropriate averaging. 

The latter consideration suggests that nuclear reso- 
nance can be used as an independent and accurate mag- 
netometer to study the long term performance of the ac- 
tive compensation system in terms of Allan variance of 
the nuclear precession frequency. This procedure requires 
the v n evaluated to be determined with a very small mar- 
gin of uncertainty 5v n . In Sec. lIVI we report a practical 
example, where the best fit procedure gives a relative 
frequency uncertainty of bv n jv n rj 1.6 x 10 -6 , meaning 
that the relative field drifts in this amount, occurring in 
the averaging time (2800 sec), are appreciable. The rela- 
tively poor signal-to-noise ratio of our setup makes this 
procedure difficult. Nevertheless it can be proposed as a 
characterization method for evaluating performance over 
the very long term, especially in conjunction with more 
sensitive magnetometers. 



C. Residual noise due to field inhomogeneities 

As discussed above, our system allows for complete 
control of all the dc components of the magnetic field, 
while tracking only its modulus over time, and counter- 
acting its variations through adjustments of the strongest 
(z) component. It is worth considering the limitations 
deriving from such a configuration, with the specific aim 
of evaluating the critical features that occur when the 
system is used in a dual channel setup. 

Let us assume that the system perfectly compensates 
the field fluctuations in the main arm, and let us derive 
the residual noise in the secondary arm. We will assume 
that the noise B n and the compensation field B c z are ho- 
mogeneous, while the bias field may differ due to static 
gradient: the field is Bq in the main cell, and B\ in the 
secondary one, with B\ w Bq + J Ax, J being the Jaco- 
bian of the field and Aa; the position of the secondary 
cell with respect to the main one. Thus, perfect compen- 



sation in the main cell reads 

(B + B n + B cZ y = B 2 (4) 
while, concerning the secondary cell, 

(Bo + Bn + B c z + JAxy =B\ + A(Bl), (5) 

A(B 2 ) representing residual noise. The difference be- 
tween eqs.Q]and[5]gives 

2(B n + B c z) -{3AS) = A(B 2 ). (6) 

Solving cq. [4] for B c in the hypothesis Bo — Bqz gives 

/ B 2 \ B 2 

B c = -B -B n -z + B ll 2i ^ 

° V B 2 J 2B 

where B n ± is the projection of the noise field on the xy 
plane. Thus, for small noise, in the first order approxima- 
tion, the first parenthesis in eq.[6] reduces to 2B± ■ 3 Ax. 
In conclusion, using a scalar sensor to compensate homo- 
geneous field noise in a dual channel setup requires accu- 
rate preliminary compensation of the gradient of the field 
components perpendicular to the bias field, in our case 
dB x /dx and dB y /dx. The inhomogeneity of the parallel 
component dB z /dx has a direct (i.e. first-order) effect 
on the field modulus (and thus must be compensated to 
force the secondary arm at the center of its resonance) 
but has no first-order effects on the noise compensation 
efficiency. 

IV. PERFORMANCE IN APPLICATION 

The compensation system was used in an ultra-low- 
field NMR experiment, with measurements similar to 
those reported in [22j |. Nuclear precession of remotely 
polarized water protons was detected in a fiT field. 
Fig. [6] shows the free induction decay (FID) signals and 
their power spectral density, obtained after averaging 700 
traces. 

Active compensation of the bias field renders unneces- 
sary the frequency rescaling procedure described in [22j], 
where NMR traces were interpolated and time-rescaled 
to compensate the nuclear Larmor frequency drift, which 
was inferred from the trace-by-trace measured bias mag- 
netic field. 

Furthermore, operating in locked-field conditions im- 
proves the quality of the average trace, as can easily be 
recognized by comparing the residual traces obtained as 
a difference of the average signals with their best fits 
(modeled as an exponentially decaying sinusoid). In the 
best-fit examples shown in Fig. [HI activating the field sta- 
bilization system leads to a reduction of the minimum \ 2 
by a factor of 7. In addition, locking the field significantly 
reduces the low frequency noise, as can be clearly seen in 
the plot (d). 
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FIG. 6: The plots on the left show FID signals (dots), fitting 
curves (thin grey lines), and residuals (thick lines), in the time 
domain. The corresponding spectra are reported on the right. 
The upper plots are obtained with active field stabilization, 
while the lower plots are obtained by tracking the field and 
applying the methods described in [22(]. Time-domain traces 
were improved (as in [23]) by identification and subtraction 
of components at 50 Hz and its multiples, and linear filtering 
(1 st order Bessel filter, 170-^190 Hz bandpass). 

V. CONCLUSION 

We have presented and characterized a simple, eco- 
nomic and reliable magnetic field compensation system. 
The system is integrated in a set-up working at a bias 
field of a few /iT, which is measured by means of a dual- 
channel magnetometer, based on all-optical atomic scalar 
sensors and operating as a self-oscillator. 

The compensation system is driven by one of the two 



channels and acts on the field component parallel to the 
bias, feeding back the current supplied to a Hclmholtz 
pair. It counteracts the common mode magnetic field 
fluctuations, thus facilitating difference-mode detection. 
The difference- mode signal reproduces the field anoma- 
lies generated by a weak source located in the proximity 
of one sensor. 

The operating principle makes the compensation par- 
tial, meaning that only the modulus of the magnetic field 
is kept constant. The limitations posed by this feature 
have been discussed, arriving at the conclusion that the 
compensation efficiency is mainly degraded by inhomo- 
geneities of the field components perpendicular to the 
bias direction, which therefore require careful control. 

Detailed data analyses were carried out to characterize 
the short-term performance of the stabilization system 
and to evaluate the phase noise spectra, from which the 
stray field extinction ratio was inferred. Concerning the 
power-net disturbances, we demonstrated an extinction 
factor in excess of 40 dB. 

As a practical example, we reported results obtained 
with and without active compensation in an ultra-low- 
field NMR experiment. The free induction decay of re- 
motely polarized water samples was recorded and aver- 
aged over long intervals. The possibility of using these 
NMR measurements to assess the long term performance 
of the compensation system was also considered. 
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